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Demodulation for Optimum Correlat ion of GPSSignals 

* 

Field of the Invention 

5 The present invention relates to a method and system for increasing the 

sensitivity of the demodulation of direct-sequence spread spectrum modulated- 
signals and more particularly to the application of a Viterbi demodulation type 
algorithm to the output of a coherent correlation of a GPS C/A LI signal to improve 
the demodulation sensitivity. 

10 

Background to the Invention 

Due to Federal Communications Commission (FCC) requirements for 
identifying mobile phone location, vendors are beginning to incorporate global 
positioning system (GPS) receivers into their handsets. This is an assisted form of 

15 GPS where the host base station provides key search parameter information to the 
handset such that the GPS satellite search and subsequent delay measurement 
becomes significantly less computationally extensive for the handset. In the code- 
division multiple access (CDMA) context, the messaging between the base station 
and the handset to support the assisted GPS (AGPS) is documented in "Position 

20 Determination Service Standards for Dual Mode Spread Spectrum Systems," 
TIA/EIA/IS-801-T, published by the Telecommunications Industry Association (TIA). 

In order for mobile location to be successfully based on AGPS, it is necessary 
that the GPS receiver in the handset unit have significantly higher sensitivity than is 
nominally required in a typical stand alone GPS unit. The reason is that the handset 

25 unit is usually oriented in positions or in locations that are not favorable for GPS 
signal reception. Hence a signal to noise ratio (SNR) sensitivity target is typically 
around 17 dB-Hz. Given that the GPS Coarse/Acquisition (C/A) L1 band which is 
typically used for mobile position location has a bandwidth of about 2 MHz, this 
implies that the input raw GPS, prior to processing, is more than 40 dB below the 

30 thermal noise floor. 

There is also a significant push to drive the cost of the overall handset down 
which limits the quantity of the processing which can be done on the received GPS 
signals to extract them from the noise. This also limits the quality of the various 
radio frequency (RF) oscillator components used in the basket's receiver. 

35 Consequently, it can be expected that there will be some additional local oscillator 
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(LO) noise that needs to be compensated for. Notably, there is no closed automatic frequency 
correction (AFC) loop possible during the GPS measurement as the input SNR is too low. The 
lack of AFC, and the use of low cost, low power, RF LO components, implies that there will 
inevitably be a significant finite offset, drift, and some instability associated with the LO down 
conversion frequency. 

The C/A signal from a GPS SV (space vehicle) is a DS-SS (direct sequence spread 
spectrum) with a chip rate of 1 .023 Mcps. It is modulated as binary phase shift keying (BPSK) 
on a 1.574 GHz carrier. The GPS receiver correlates the received signal with a locally 
generated DS-SS code signal, f n the AGPS scheme, traditional DS-SS correlation is also done. 
However, detailed information regarding the doppler shift of the SV GPS signal and code offset 
is available from the host base station (BS) which significantly reduces the search and detection 
effort. Nevertheless, the mobile GPS receiver is still required to determine the code delay and 
doppler to a finer resolution than that available from the host BS such that mobile location is 
possible. 

Based on standard assumptions regarding the noise in the GPS signal channel, the 
optimum receiver would correlate the signal in a coherent fashion over an integration time 
period that is sufficiently long to provide about 11 dB SNR at the correlator output This will 
typically provide an adequate probability of detection with a reasonable false alarm rate. 
However, in the case of a mobile GPS receiver, due to the instability of the RF LO and the 
uncertainly of the SV doppler, the coherent integration epoch needs to be limited. Also valid 
GPS readings are still required even if the user does not hold the receiver steady. Hence, . 
typically the coherent integration time is limited to 1 0 msec or less. As the available coherent 
integration epoch is not sufficient to obtain the sensitivity required, non-coherent summations of 
sequential coherent correlation integration outputs are used. However, non-coherent 
processing is a very inefficient means of further enhancing the SNR of a signal as rt discards 
certain Known statistical aspects of the signal. 

Various examples of the use of Viterbi algorithms for phase trajectory determinations are 
known. For example US Patent No. 6.477,208 to Huff teaches a method and apparatus for 
processing a received digitally-modulated carrier signal to coherently demodulate a signal 
utilizing a composite trellis diagram. However, the method of Huff is only applicable to cases of 
discrete jumps in phase based on specific modulations. When using a continuous, free running 
oscillator, the phase steps are continuous in time. Huff does not teach, a method for 
approximating these phase steps. 
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ST Summary of the Invention 

The method and system of the present invention overcome the above by providing a 
way to demodulate the RF phase trajectory of a DS-SS correlation, using coherent integration 
for a small time epoch. The receiver fn the present invention downshifts an incoming signal 

5 using a local oscillator. The receiver then despreads the signal. 



Empf.zeit: 19/10/2004 19:16 



2a 

Empf .nr.: 807 P. 007 



WO 2004/049663 PCT/CA2002/001843 

In a perfect system the output from the despreader would be a DC signal. 
However, due to frequency and phase instabilities in the radio frequency local 
oscillator, and due to uncertainty of the GPS SV doppler, the phase trajectory is 
generally a time varying phasor. 
5 In order to provide increased sensitivity, a trellis is created with a fixed 

number of phase states evenly distributed between zero and 360°, The state 
transitions in the trellis are limited based on known phase slew rate limitations. 
Probabilities for these state transitions are similarly based on these same factors. 

A Viterbi algorithm is then applied to the phase of the despread signal. 

10 Likelihood metrics are created for each path based on the phase of the despread 
signal as well as path probabilities. The system and method then uses the Viterbi 
algorithm to find the highest likelihood path for the phase irajectory through the trellis. 

The present invention therefore provides a method of improving sensitivity in 
the demodulation of a received signal over an . arbitrary measurement time epoch, 

15 said method comprising the steps of: correlating said received signal in a coherent 
fashion; and utilizing a Viterbi phase state keying trellis demodulation with a variable 
resolution of phase states over 0 to 360° to demodulate the radio frequency phase 
trajectory of said signal throughout the measurement time epoch. 

The present invention further provides a method of increasing sensitivity of a 

20 receiver for direct-sequence spread spectrum signals, said method comprising the 
steps of: downconverting the direct-sequence spread spectrum signal; despreading 
the downconverted signal creating a despread signal with a radio frequency phase 
trajectory; and applying a Viterbi algorithm to the despread signal to determine the 
radio frequency phase trajectory , whereby the algorithm has the steps of: defining a 

25 number of possible states, each state corresponding with a discrete phase and all the 
states together creating an even distribution of phases between zero and 360°; 
creating possible state transitions from each state; setting the probability of a given 
state transition occurring; defining a trellis with paths based on the possible state 
transitions at a number of discrete time intervals; taking a sample of the despread 

30 signal at a number of discrete time intervals; creating a likelihood metric for each 
path based on the phase trajectory of the despread signal at the time interval of the 
path; discarding all paths entering a state except for the path with the highest 
likelihood metric entering the state; and determining the most likely radio frequency 
phase trajectory based on the paths not discarded. 

35 The present invention still further provides a receiver for receiving a direct 

sequence spread spectrum signal, said receiver comprising: an antenna for receiving 
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the direct sequence spread spectrum signal; a downconverter for downconverting the 
received signal, producing a downconverted signal; an analog to digital converter to 
convert the downconverted signal to a digital signal; a despreader for despreading 
and coherently correlating the digital signal to a known signal, creating a despread 
5 signal; and a processor for applying a Viterbi algorithm to the despread signal. 

Brief Description of the Drawings 

Preferred embodiments of the present invention will now be described in 
greater detail and will be better understood when read in conjunction with the 
10 following drawings in which: 

Figure 1 is a schematic diagram of the receiver of the present invention; 

Figure 2 is an example trellis diagram for use with the method of the present 
invention; 

Figure 3 is an alternative trellis* diagram for use with the method of the 
15 present invention; 

Figure 4 is a graph illustrating an approximation of the probability distribution 
function of a phase sample as a periodic gaussian pulse on top of a constant 
function; and 

Figure 5 is a graph illustrating the losses in decibels verses the signal to noise 
20 ratio of the samples received in decibels. 

Detailed Description of the Drawings 

In this disclosure, a system and method are described that provide a practical 
way of performing an optimal correlation of the DS-SS GPS signal using the known 

25 statistical properties of the signal. In particular, the method makes use of the 
generally known slew rate limitation of the LO instability and SV doppler. 

The main aspect of the method is a Maximum Likelihood Viterbi procedure 
similar to the processing used in trellis decoding of standard PSK communication 
signals. Although the method is applied exclusively to GPS signals in the following 

30 disclosure, one skilled in the art will realize that the method is also applicable for 
precision tracking of CDMA pilot signals for applications of high sensitivity mobile 
position location. 

Reference is now made to Figure 1 . The antenna 10 of a receiver 1 receives 
a direct sequence spread spectrum signal. The signal is then passed through band- 
35 pass filter 20 to isolate the signal. In a preferred embodiment, band-pass filter 20 is 
centred at 1 .574 Ghz, which is the frequency of the GPS C/A L1 carrier. 
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The output from band-pass filter 20 is then downconverted by multiplying it 
with the signal from the radio frequency local oscillator 30 in multiplier 40, producing 
a baseband signal. This baseband signal is then passed through low-pass filter 50 
and an analog to digital converter 60 to create a spread spectrum digitized signal. 
5 The spread spectrum digitized signal is then passed through a digital signal 

processor 70. Digital signal processor 70 is comprised of a despreader 80 and a 
processor 90. Despreader 80 is comprised of multiplier 84 and signal generator 83 
which generates coarse acquisition LI signal 82 which is multiplied with the spread 
spectrum digitized signal. 

10 In a perfect system, the output from despreader 80 would be a DC signal. 

However, due to frequency and phase instabilities in local oscillator 30, and due to 
uncertainty of the GPS SV doppler, the phase trajectory is generally a time varying 
phasor. The output of despreader 80 is thus processed by processor 90, which 
applies a Viterbi algorithm as described below. 

15 Consider the total GPS signal measurement time epoch, denoted by T to t as 

being divided into K equal length segments of T c each. It is assumed that the GPS 
correlator will coherently integrate for each of the segment lengths T c resulting in K 
correlation samples X=[x 1> x 2 ,...,x K ] T . It is also assumed that T c is short enough such 
that uncertainties in GPS doppler, RF LO offset or frequency/phase instabilities of 

20 the RFLO are such that the resulting signal loss in the coherent correlation output 
samples in X is not significant. Typically T c will be somewhere between 1 to 10 
msec. 

X=[x 1l x 2l ...,x K ] T is mapped into corresponding sample phases denoted by 

P=tPi,P2 Pk] t - These phases will be uniformly distributed between 0 and 360°. In 

25 order to implement a practical Viterbi algorithm, discrete phases are defined as being 
uniformly spaced such that: 

M = number of discrete phases 

So,s 1 ,...s M -i are the discrete phase states given by s n =(m/M) 360° 
Let k denote the k m state interval corresponding to the k m sample time or x k . 
30 Reference is now made to Figure 2, which shows an example trellis diagram for 
better illustration of the following description. 

Consider the m m state of the k m interval. If there are no phase dynamics due 
to the RF LO or doppler, then it will map into state m of the (k+1) lh interval with a 
probability of 1. However, due to the uncertainty in the phase dynamics, the m m state 
35 can map into state (m-1), m, (m+1) of the (k+1) th interval. Clearly, with these trellis 
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state transition probabilities, the modeled slew rate of the phase dynamics can be 
accurately limited. 

The trellis of Figure 2 is only one example. One skilled in the art will realize 
that other trellis state transition probabilities are possible, and that a probability trellis 
5 can be created to more accurately model these other probabilities. For example, in 
Figure 3, a system is illustrated in which the slew rate can be higher but a lower 
probability is assigned to higher slew rates. 

Intuitively such slew rate limiting lattice structures as described above 
effectively "low pass filter" the phase samples to reduce the influence of noise in 
10 correlation samples X. 

The number of phase states, "M", is a function of the phase resolution 
required. If M is too small then there will be a small demodulation loss due to the 
discrete phase quantization. 

At each update, the standard Viterbi algorithm is used to eliminate the trellis 

15 paths of lower likelihood. (A good description of the Viterbi algorithm is found in 
"Digital Communications" by J. Proakis, McGraw Hill , 3 rd edition 1995.) Briefly, at 
the k lh update a finite number of paths converge onto each of the M trellis phase 
states. For example, in the trellis structure of Figure 2, there are three candidate 
paths and in the trellis structure of Figure 3, there are five such paths. Each of these 

20 paths will have a likelihood metric associated with it. The path selected for each path 
is the one associated with the highest likelihood. Consequently, at each interval 
step, P paths converge to each trellis state and subsequently P new path segments 
are generated, as shown in the above examples for the (k+1) lh trellis interval. 

The measured samples, p k , are the phases of the complex correlator outputs 

25 x k =x ik + iXq k as determined by: 

p k = atan(x qk /x ik ). 

This is an example only as there are a variety of numerical methods of 
extracting the phase from x k . In order to update the likelihood metric for each 
surviving path at the k m trellis interval, an estimate of probability of the received 

30 phase being p k conditioned on the phase state and the trellis path corresponding to 
the particular candidate path being considered is required. As the SNR of the 
correlator outputs, x kl is typically low (usually less than 0 dB), the probability 
distribution function (PDF) of the phase sample p k is not gaussian. Also the phase 
trajectory is dependent in a complicated statistical way on the candidate trellis path 

35 being considered. Consequently, finding the actual conditional probability is 
impractical. However, there are several assumptions and approximations that can be 
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reasonably made that significantly simplify the likelihood probability metric. These 
are: 

1 . Model the phase random process as first order Markov 

2. Approximate the PDF of p k as a periodic gaussian pulse on top of a 
5 constant function as illustrated in Figure 4. 

Given this approach, the conditional probability metric can be approximated 

as (Pk-Sm) 2 where S m represents the assigned phase of the m th state. The influence 

of the constant portion of P(Pk|S m ) is irrelevant as it is added to each candidate trellis 

path and therefore does not result in a path differentiating factor. Practical simulation 
» 

10 experience has shown that this simple path likelihood metric is effective even though 
the approximation to the conditional probability of P(pk|S m ) is crude. 

Another very important factor is that the differentiating part of the likelihood 
metric when the above approximation is used becomes independent of the variance 
of the conditional probability of p k . This is a key factor in the robustness and 

15 practicality of the method. 

« 

Example 

Below there is provided a simple MatLab simulation to demonstrate the 
effectiveness of the present method. 

20 In the simulation, an array of 2000 complex samples of x k were generated 

based on the assumption that the signal component was a constant amplitude 
exponential with a phase increment of 0.03 radians. Random noise was added to 
the samples of x k that was zero mean gaussian to both in-phase and quadrature 
phase components. The noise components were independent from sample to 

25 sample and also the noise in the quadrature and in-phase components were 
independent of each other. The variance of the noise components was selected to 
satisfy a given SNR of the correlation samples x n for the simulation. 

The Viterbi demodulator was implemented as described above and tracked 
the phase excursion of the correlation samples x k . The accuracy of the demodulated 

30 phase was compared to the actual phase to determine a signal loss factor that was 
dependent on phase errors only. 

Figure 5 shows a plot of this loss factor as a function of the SNR of the 
correlation samples. As expected, the loss generally decreases as a function of 
increased SNR. At very low SNR levels, the phase tracking becomes so erratic that 

35 the Viterbi demodulation is not effective. . 
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Figure 5 can be used to determine the number of samples required for a 
given SNR. Consider the samples at 0 dB SNR. As stated earlier, to satisfy the 
requirements of adequate detection and false alarm statistics, approximately 11 dB of 
SNR is required at the output of the overall correlation process. Hence if the 
5 samples Xr of 0 dB SNR are generated after coherent correlation of T c , the number of 
samples required, K, is given by the formula: 

10log(K)>= (Minimum Correlation Process Output SNR) - (Coherent 

Correlation SNR) + (coherence integration losses in dB) 

or, using Figure 5, 
10 10log(K) >= 1 1 - (0) + 1 = 12 dB. 

Solving this, K=1 6. 

As indicated, the additional 1 dB added to the right hand side of the above 
formula is to offset the coherence integration losses represented in Figure 5. 
Comparing this with the traditional scheme which uses non-coherent summation of 
15 the correlator output samples, a minimum of K=26 is required. Hence, for a given 
GPS hypothesis dwell time, sensitivity can be significantly better with the Viterbi style 
demodulator as represented in this disclosure. 



20 Matlab simulation listing 

% viterbi demodulator for phase excursion 
% 

25 % number of state transitions 
K = 2000; 

% number of states 
M = 30; 

30 

% phase angle of each state 
state_angle = (2*D/M) * [Q:M-1]'; 

% model the actual signal phase 
35 phi_sig = .03 * [1:K]'; 

% ac t ua 1 s igna 1 
SNRdB = -6; 
SNR = 10^(0.1* SNRdB ) ; 
40 u = exp (1 j*phi_sig) ; 

x = randn(K,l) + 1 j *randn (K, 1) + sqrt(SNR) * u; 
phi_meas = angle (x) ; 

■ 

% metric for phase state 
45 % at each update there is only one surviving path for each 
state 
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Mt = zeros (M, 1) ; 

% paths, and initiate the k=l state; equal probability that we 
can start from 
5 % any state 

P = zeros <M, K) ; 
P( : , 1) = [1:M] 1 ; 
Pnew = P; 

10 % loop over the overall set of K state transitions 
I = ones (M, 1) ; 

for k=2:K 

15 % determine the path metric of the new segment from state 

k-1 to state k 

gl = phi — meas(k) * ones(M,l) - state_angle; 

g2 = gl + 2*D*I; 

g3 = gl - 2 * D*I; 
20 M_new = (min(abs( [gl,g2,g3] ' ) ) ' ) ."2; 

% determine the new path for state m 
for m=l:M 

if m == 1 
25 mml = M; 

m0 = 1; 
mpl = 2; 
elseif m == M 
mml = M-l ; 
30 m0 = M; 

mpl = 1 ; 

else 

mml = m-l; 
mO = m; 
35 mpl = m+l ; 

end 

[gl,i] = min( [Mt (mml) ,Mt (mO) ,Mt (mpl) ] ) ; 
if i == 1 

40 Pnew(m # l:k) = [P (mml , 1 :k-l) ,m] ; 

Mt(m) = Mt(mml) + M_new(m); 
elseif i ==2 

Pnew(m,l:k) = [P(m0, l:k-l) ,m] ; 
Mt(m) = Mt(m0) + M_new(m) ; 

45 else 

Pnew(m,l:k) = [Ptmpl, 1 :k-l) ,m] ; ■ 
Mt(m) = Mt(mpl) + M_new(m) ; 

end 

end 

50 

P = Pnew; 



55 end % end of the K transition periods 
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figure ( 1 ) ; plot ( P ' ) ; 
hold on; 

[gl, i] = min(Mt) ; 
5 plot(P(i, : ) , 'or' ) ; 
hold off; 

% unwrap P 
Puw = zeros (K, 1) ; 
10 Puw(l) = P(i, 1) ; 
offset = 0; 
for k=2:K 

Puw(k) = P(i,k) + offset; 
if Puw(k) - Puw(k-l) > 1 
15 offset = offset - M; 

Puw (k) = Puw(k) - M; 

end 

if Puw(k) - Puw(k-l) < -1 
offset = offset + M; 
20. Puw(k) = Puw(k) + M; 

end 

if abs(Puw(k) - Puw(k-l)) > 1 

error ( 1 discontinuity in Puw'); 

25 end 
end 



figure (2) ; 
30 plot (Puw); 

figure (3 ) ; 
plot(P(i, :) , 'r' ) ; 
hold on; 

35 

plot (mod (phi_sig * M/(2 * □),M), , g'); 
hold off; 
grid on; 

40 

% calculate the integration loss 

Loss = 20 * loglOlabslu 1 * exp(lj*state_angle(P(i, : ) ) ) ) /K) 



45 

The above-described embodiments of the present invention are meant to be 
illustrative of preferred embodiments and are not intended to limit the scope of the 
present invention. Also, various modifications, which would be readily apparent to 
one skilled in the art, are intended to be within the scope of the present invention. 
50 The only limitations to the scope of the present invention are set forth in the following 
claims appended hereto. 
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